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1 Introduction 

The inverse problem of Lagrangian mechanics is the question: given a system of second- 
order ordinary differential equations, under what circumstances does there exist a regular 
Lagrangian function, such that the corresponding Lagrange equations are equivalent (i.e. 
have the same solutions) as the original equations. Locally, the question can be translated 
immediately into more precise terms as follows: considering a given second-order system 
in normal form 

q' = r{q,q), (1) 

which (for the time being) we take to be autonomous for simplicity, what are the con- 
ditions for the existence of a symmetric, non-singular multiplier matrix gij{q,q) such 
that 
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for some L. Clearly {gij), if it exists, will become the Hessian of the Lagrangian L. 
The literature on this problem is extensive; the conditions for the existence of L are 
usually referred to as the Helmholtz conditions, but these can take many different forms 
depending on the mathematical tools one uses and on the feature one focusses on. For 
a non-exhaustive list of different approaches see [24], [5], [26], [22], [2], [9], [15], [1], 
[19], [3]. In this paper, the tools stem from differential geometry and therefore provide 
coordinate-free results. In addition, while we will actually study generalizations of the 
above problem which allow for certain classes of non-conservative forces, the attention 
will be mainly on conditions on the multiplier g. 

We will consider two types of non-conservative forces, leading to Lagrangian equations 
of one of the following forms: 

dt \dq^ ) dq^ dq^' ^ ' 

or 

d f dL\ dL t 

d-tWO'dt^"""'^'^'' ^fc^ = -^^^- (3) 

In the first case, when the function D is quadratic in the velocities (and —D is positive 
definite) the classical terminology is that we have dissipation of Rayleigh type (see e.g. 
[13]); we will not put restrictions on the form of D, however. In the second case, in which 
the existence of a D as in (2) is excluded, the right-hand side is often referred to as a 
gyroscopic force (sec e.g. [23]). 

Perhaps we should specify first what we will not do in this paper. In older contributions 
to the inverse problem for dissipative systems, such as [10], the emphasis was on trying to 
recast a dissipative system into the form of genuine Euler-Lagrange equations, that is to 
say that in the case of given equations of type (2) one would try to find a different function 
L' such that the Euler-Lagrange equations of L' are equivalent to the given system. In 
contrast, our goal here is to study under what circumstances a given second-order system 
in normal form (1) can be recast into the form (2) (or (3)) for some functions L and D 
(or L and ujki)- 

In order to explain our objectives in more precise terms, let us recall first some of the 
different ways of characterizing the inverse problem conditions in the classical situa- 
tion. The natural environment for a second-order system is a tangent bundle TQ, with 
coordinates {q, v) say, where it is represented by a vector field F of the form 

r = »'!- + /'(,,») A. (4, 

If 5 = {d/dv"^) ^dq^ denotes the type (1, 1) tensor field which characterizes the canonical 
almost tangent structure on TQ [6, 14], F represents a Lagrangian system provided there 
exists a regular Lagrangian function L such that (see e.g. [27]) 

CT{S{dL)) = dL- (5) 
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9l ■= S{dL) is the Poincare-Cartan 1-form. The above condition is perhaps the most 
compact formulation of the problem, but has little or no practical value when it comes to 
testing whether such an L exists for a given F. A shift of attention towards the existence 
of a multiplier leads to the following necessary and sufficient conditions [5] : the existence 
of a non-degenerate 2-form u G /\^{TQ), such that 

Croo = 0, tj(X^,F^) = 0, = 0, VX, F, Z G A'(M). (6) 

Here and X" refer to the vertical and horizontal lift of vector fields, respectively. 
The latter makes use of the canonical Ehresmann connection on r : TQ — ^ Q associated 
with a given second-order vector field T: in coordinates, the vertical and horizontal lift 
are determined by 



_dj^_d_ H .= — = — -V^— hr r^=-i^ 
dq^ dv'' ' ■ dq' dq' ' dv^ , w ere • ^ , 



Such a 2-form oj will be closed, hence locally exact, and as such will be the exterior 
derivative dOi for some Lagrangian L. At this point it is interesting to observe that the 
2n X 2n skew-symmetric component matrix of lu is completely determined by the n x n 
symmetric matrix 

_ d^L 

The matrix (gij) geometrically represents the components of a (0,2) symmetric tensor 
field g along the tangent bundle projection r, and the relationship between u on TQ and 
g along r has an intrinsic meaning as well: oj is the Kahler lift of g (see [21]). A more 
concise formulation of the Helmholtz conditions therefore, when viewed as conditions on 
the multiplier g, makes use of the calculus of derivations of forms along r, as developed 
in [20, 21]. We will show in the next section how both the conditions (5) and (6) have 
an equivalent formulation in those terms, and this will be the basis for the generalization 
to Lagrangian systems with non-conservative forces, which will be the subject of the 
subsequent sections. 

The first authors to discuss the inverse problem, in the sense of analyzing the conditions 
which a given representation of a second-order system must satisfy to be of the form 
(2), were Kielau and Maisser [18]. We showed in [8] how the results they obtained 
via an entirely analytical approach can in fact be reduced to a smaller set. But we 
also argued in the concluding remarks of that paper that the more important issue 
is the one we formulated above, which starts from a normal form representation of the 
dynamical system. For that purpose it is better to approach the problem in a coordinate- 
independent way, i.e. to make use of the tools of differential geometry already referred to. 
We will see that the methods we will develop for the dissipative case (2) apply equally 
to the gyroscopic case (3). To the best of our knowledge the latter problem has not been 
dealt with before in its entirety (though a relevant partial result has been published in 
[19]). An additional advantage of the coordinate- independence of our conditions is that 
they cover without extra effort results such as those derived in [17] for the description 
of Lagrangian systems in 'nonholonomic velocities'. In Section 3 we follow the lines of 
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the construction of Helmholtz conditions on the multipher g for the standard inverse 
problem, and arrive in this way at necessary and sufficient conditions which involve g 
and D in the dissipative case, and g and lo in the gyroscopic situation. At the end 
of this section we briefly discuss how the partial result mentioned above is related to 
our work. In Section 4 we succeed in eliminating the unknown D and 00 altogether 
to arrive at necessary and sufficient conditions involving the multiplier g only. This is 
particularly interesting, because a given F may actually admit multiple representations 
of the form (2) for example. In other words, different choices of a multiplier g may exist, 
which each require an adapted (generalized) dissipation function D to match the required 
format. In fact it cannot be excluded that a given F may actually have representations 
in the form (2) and (3) at the same time, of course with different multipliers g (and 
thus different Lagrangians L). We will encounter such situations among the illustrative 
examples discussed in Section 5, where we also briefly indicate in the concluding remarks 
how the whole analysis can be carried over to the case of time-dependent systems. In 
an appendix we make an excursion to a different geometrical approach which in fact is 
essentially time-dependent: we use techniques from the theory of variational sequences 
to relate our results more closely, at least in the dissipative case, to those obtained in 
[18], which after all was the work which first brought this subject to our attention. 



2 Basic set-up 

In order to keep our analysis reasonably self-contained, we need to recall the basics of 
the calculus of derivations of forms along the tangent bundle projection r : TQ Q. 
Vector fields along r are sections of the pull-back bundle t*TQ — ^ TQ and constitute a 
module over C°°{TQ)^ denoted by X{t). Likewise, a fc-form along r assigns to every point 
Vq € TQ an exterior /c-form q = T{vq) G Q; wc use the symbol /\(t) for the C°°{TQ)- 
module of scalar forms along r and V{t) for the module of vector- valued forms. The 
theory of derivations of such forms, as established in [20, 21], follows closely the pioneering 
work of Frolicher and Nijenhuis [11]. The difference is that there is a natural vertical 
exterior derivative available, but a full classification requires an additional horizontal 
exterior derivative d" , which must come from a given connection: in our situation, this 
is the connection associated with F mentioned earlier. We limit ourselves here to a brief 
survey of the concepts and properties we will need. An elaborate version of the theory 
(with rather different notations) can also be found in [30]. 

Elements of /\ (r) in coordinates look like forms on the base manifold Q with coefficients 
which are functions on TQ. Thus they can be seen also as so-called semi-basic forms 
on TQ, and we will generally make no notational distinction between the two possible 
interpretations. It is clear that derivations of such forms are completely determined by 
their action on C°°{TQ) and on /\^{Q)- As such, the vertical and horizontal exterior 
derivatives are determined by 

dT =Vi{F)dq\ d"F = Hi{F)dq\ FeC^{TQ), 
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(Tdq' = 0, d"dq' = 0. 

Obviously, for L G C'^{TQ), d^L e /\ (t) has the same coordinate representation as 

S{dL) € /\^{TQ); in Hne with the above remark therefore, we will also write 9^ = d^ L 
for the Poincare-Cartan 1-form. Derivations of type are defined as in the standard 
theory. For A ^V{t), we put 

dl = [iA,d^], d'i = [i^,d% 

and call these derivations of type d"^ and respectively. The action of all such deriva- 
tions can be extended to vector-valued forms and then another algebraic type derivation 
is needed for a classification, but we will introduce such extensions, which can all be 
found in [20, 21], only when needed. The horizontal and vertical lift operations, already 
referred to in the introduction, trivially extend to vector fields along r and then every 
vector field on TQ has a unique decomposition into a sum of the form X" + , with 
X,Y ^ X{t). Looking in particular at the decomposition of the commutator 
suffices to discover two important derivations of degree zero: 

They extend to forms by duality and are called the horizontal and vertical covariant 
derivatives. In coordinates 

B^^F = X'Vi{F), Dx^ = 0, D^dg^ = 0, 

Df F = X^H,{F), Df A = x%{r^)-^, Df dg^ = -X^V,{rl)dq\ 

For later use, we mention the following formulas for computing exterior derivatives of, 
for example, a 1-form a or a 2-form p along r: 

d^a {X, Y) = D^a (Y) - D^a (X), a G A'W, (8) 
d^p {X, Y,Z)=J2 ^xP Z), p G A'(r), (9) 

X,Y,Z 

and similarly for d". Here J2x YZ represents the cyclic sum over the indicated arguments. 
It is also of interest to list the decomposition of the other brackets of lifted vector fields: 

[X^,Y^] = (D^Y-D^X)'', 

[X", y"] = (Df y - B^Xf + {R{X, Y)y . 

The latter relation is just one of many equivalent ways in which the curvature tensor 
R e F^(r) of the non-linear connection can be defined. The connection coming from F 

has no torsion (since (7) obviously implies that Vi{Tj^) = V^.(r^)): it follows that d^ and 
d" commute. In fact the commutation table of the exterior derivatives, for their action 
on scalar forms, is given by 

=d^d^ = 0, d^d« = -d^d^, ^[d'',d"]=d"d'' = d]i. (10) 
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Finally, the given dynamical system T comes canonically equipped with two other oper- 
ators which are crucial for our analysis, namely the dynamical covariant derivative V, a 
degree zero derivation, and the Jacobi endomorphism $ G V^{t). Again, the simplest 
way of introducing them comes from the decomposition of a Lie bracket: they are the 
uniquely determined operations for which, for each X G X{t), 

[r,x"] = (vx)" + {^xY. 

The usual duality rule V(X, a) = (VX, a) + (X, Va) is used to extend the action of V 
to 1-forms, and subsequently to arbitrary tensor fields along r. In coordinates, 

V^ = r(n v(|-)=rjA, V{dq^) = -T]dq\ (11) 

and 

^ = -^-^'n-m). (12) 

One clear indication of the importance of these operators is the following link with the 
curvature of the connection: 

d^# = 3i2, d"<^ = VR. (13) 



We are now ready to go back to the generalities about the inverse problem discussed in 
the previous section. To begin with, using the tools which have just been established, the 
compact formulation (5) of the inverse problem is equivalent (see [21]) to the existence 
of a regular function L G C°°{TQ) such that 

VOl = d'^L. (14) 

Secondly, the necessary and sufficient conditions (6) now really become conditions on 
the multiplier matrix; they are translated via the Kahler lift [21] into the existence of a 
non-degenerate, symmetric (0, 2)-tensor g along r satisfying the requirements 

Vg = 0, g{^X,Y) = g{X,^Y), B^xg{Y, Z) = B^g{X, Z). (15) 

It is possible to prove directly that (14) implies (15) and vice versa (a sketch of such a 
proof was presented in [4]). We will not show how to do this here, however, as it can 
easily be seen later on as a particular case of the more general inverse problem studies 
we will start analyzing now. 



3 Lagrangian systems with dissipative or gyroscopic forces 

Consider first equations of type (2). It is obvious that, at the level of a characterization 
like (5) , a given second-order field F will correspond to equations of type (2) if and only 
if there exist a regular function L and a function D such that 

£riS{dL)) =dL + S{dD). (16) 
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We take the opportunity to illustrate first how such a relation, when stripped to its 
bare essentials, i.e. when one observes that it is in fact a condition on only n of the 2n 
components, is transformed into a corresponding generalization of (14). To this end, 
note first that there exists a dual notion of horizontal and vertical lifts of 1-forms, from 
A^(r) to /\^{TQ), defined by a^(X«) = a{X), a^(X^) = 0, and likewise for a^. We 
then have the following decompositions, for any a G /\^{t) and L G C°°{TQ): 

Cra" = + (Va)^, (17) 
dL = {d^L)" + {d'^Lf. (18) 

Coming back to the notational remarks of the previous section: the horizontal lift is 
technically speaking the rigorous way of identifying a 1-form along r with a semi-basic 
1-form on TQ. So, when convenient, as will be the case in establishing the next result, 
we can also write 9l = (d^L)" , for example. 

Proposition 1. The second-order field T represents a dissipative system of type (2) if 
and only if there is a regular function L G C°°{TQ) and a function D G C°°{TQ) such 
that 

V9L = d"L + d''D. (19) 



Proof. We have that S{dD) = {d^ D)" for any function D, and in particular S{dL) = 
{d^ L)" = 9l. Using the decompositions (17) and (18), the condition (16) then immedi- 
ately translates into (19). □ 

Corollary 1. The condition (19) on the existence of functions L and D is equivalent to 

d"eL = 0, (20) 

which is a necessary and sufficient condition on L only. 



Proof. Using the commutator property [V, c/^] = —d" to re-express VOl = Vd^L in 
(19), we immediately get the expression 

c/"L= id^(r(L) -D), 

which is equivalent to saying that d"L = d^G for some function G. This in turn, in view 
of (10) and the triviality of d^-cohomology, is equivalent to d^Oi = —d^d^L = 0. □ 



We now want to translate these results into conditions on the multiplier g which generalize 
(15). As we observed earlier, this g will be the Hessian of L, so we look first at the relation 
between a function and its Hessian in intrinsic terms. To that end, we introduce covariant 
differentials and defined as follows: for any tensor field T along r and X G A'(r), 

D^r(x,...) = D^r(...), 
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and similarly for D^. Then for any F G C°°{TQ) we can write for the corresponding 
Poincare-Cartan 1-form Op = dXF = D^F, and define the Hessian tensor gp of F as 
gp = D^D^F, which means that 

gpiX, Y) = B^^B^F - D^y^F = B^^OpiY). (21) 

Lemma 1. For any F € C°°(TQ), its corresponding Hessian tensor gp is symmetric 
and satisfies B^gp{Y, Z) = BygpiX, Z), i.e. B^gp is symmetric in all its arguments. 
Conversely, any symmetric g along r for which B^g is symmetric is the Hessian of some 
function F. Secondly, if ^ represents any type (1,1) tensor field along t, we have 

$j gp - ($j gpf = iav^Op - d^'i^Op, (22) 

where {^a gp - {^a gpY){X,Y) := gp{^X,Y) - gp{X,^Y). 



Proof. The symmetry oi gp follows directly from d^9p = d^d^F = 0. The symmetry of 
B^gp can easily be shown by taking a further vertical covariant derivative of the defining 
relation oi gp and using the commutator property 

[B^x,'D^]=%VY-%yx- (23) 

The converse statement is obvious from the coordinate representation of the assumptions. 
Finally, making use of (8) we have 

d''i^epix,Y) = B^xiHOp)){Y) -B^mep)){x) 

= (D^$(y) - B^^(x), Op) + B\ep{^Y) - B'^ep{^x) 

= iav^ep{X, Y) + gp{X, ^Y) - gp{Y, $X), 
from which the last statement follows. □ 



We are now ready to state and prove the first main theorem, which provides the transition 
of the single condition (19) to equivalent conditions involving a multiplier 5, in precisely 
the same way as (14) relates to (15). 

Theorem 1. The second-order field T represents a dissipative system of type (2) if and 
only if there exists a function D and a symmetric type (0, 2) tensor g along r such that 
B^g is symmetric and g and D further satisfy 

Vg = D^D^D, (24) 
$j 3 - ($J gf = d^d"D, (25) 

where $ is the Jacobi endomorphism ofV. 
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Proof. Suppose T represents a system of type (2) . Then we know there exist functions L 
and D such that (19) and (20) hold true. Define g = D^D^L or cquivalcntly g{X,Y) = 
'D'^9l{Y). Obviously, g and Y)^g are symmetric by construction. Acting with V on 
and using the commutator property [V,D^] = — D^, we get for a start 

where the last line follows from (19). Now the commutator of vertical and horizontal 
covariant differentials (see [21] or [9]) is such that, at least on functions, D^D" L{X, Y) = 
B"D^L{Y,X). But 

D"D^L(y,X) - D^D^L(X,y) = d"eL{Y,X) = 0, 

in view of (20), so that (24) follows. When acting finally with on (19), we have to 
appeal to the formula for d"d^ in (10) and further need the commutator of V and d^ , 
which for the action on the module /\(t) of scalar forms is given by 

[V,d"] = 2iR + dl. (26) 

It is then straightforward to check, using (20) and the first of (13), that we get 

d^'i^Oi - idv^OL = d^d^'D, (27) 
from which (25) follows in view of the last statement in Lemma 1. 

Conversely, assume that g and D satisfy the four conditions stated in the theorem. It 
follows from the symmetry of g and g that g is a Hessian: g = D^D^F say. The 
function F of course is not unique and the idea is to take advantage of the freedom in F 
to construct an L which will have the desired properties. This is not so difficult to do by a 
coordinate analysis. Keeping the computations intrinsic is a bit more technical, but will 
give us an opportunity to recall a few more features of interest of the calculus of forms 
along r. Observe first that VD^^ is obviously symmetric, and that the same is true for 
D^Vc/ = D^D^D^D. It follows from [V, D^] = -D" that D^^? is also symmetric. Hence 

B"g{X,Y,Z) = D"D''eFiX,Y,Z) = D"D'' OpiX, Z,Y). 

If we interchange and in the last term, there is an extra term to take into account 
(since the action is on a 1-form this time, not a function). Indeed, wc have 

d«d^^f(x, z, y) = d^d^^f(z, X, Y) + eF{e{x, z)y). (28) 

Here is a type (1,3) tensor along r which is completely symmetric (and could in fact be 

defined by the above relation): its components in a coordinate basis are 0^^^ = \^V/(r^). 
Using the above two relations, expressing the symmetry of g in its first two arguments 
now leads to 

= d^d^^f(z, X, Y) - d^d^^f(z, y, X) = B^{d"eF){x, y). 
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This says that d^Op is a basic 2-form, i.e. a 2-form on the base manifold Q. On the other 
hand, we have 

where we have used the property F{X, Y) = D^D^F(F, X) again in the transition 

to the last expression. But since d"6p is basic, we can write it as U^iTd^Op, where T is 
the canonical vector field along r (the identity map on TQ), which in coordinates reads 

T = »'|,. (29) 

It follows that we can write the last relation in the form 

:= B^iVOp - d^F - iTd^Op - D^D) = 0, 

which defines another basic form /3. We next want to prove that the basic forms /3 and 
d^Op are actually closed in view of the final assumption (25) or equivalently (27), which 
has not been used so far. Keeping in mind that d" is the same as the ordinary exterior 
derivative for the action on basic forms, we easily find with the help of (10) that 

d"d"9p = yidv^Op = \d^{d^i^ep - d^d"D) = 0. 

Secondly, using also (26), 

d"l5 = d^VOp - d"d"F - d"iTd"ep - d^d^'D 

= Vd'^Op - 2iRep - dlOp - iRd^'F - df^d^Op - d^d'^D 
= Vd'^Op - idv^Op + d^'i^Op - di^d"ep - d"d''D 
= Vd^'Op - d^d"9p. 

But this is zero also because the operators V and df^ coincide when they are acting on 
basic (scalar) forms. It follows that, locally, d"6p = d"a and P = d" f, for some basic 
1-form a and basic function /. The defining relation for /? then further implies that 

d^^D = V{ep -a)- d"{F - iTa + /). 

Putting L = F — iT^a + /, the difference between L and F is an affine function of the 
velocities, so both functions have the same Hessian g, and also 9l = D^L = 9p — a. It 
now readily follows that the relation (19) holds true, which concludes our proof in view 
of Proposition 1. □ 

It is worthwhile listing the coordinate expressions for the necessary and sufficient con- 
ditions of Theorem 1. They call for a (non-singular) symmetric matrix gij(q,v) and a 
function D{q, v) such that 



VkiOzj) = V,{gik) (30) 
-gikT'; -gjkT'^ = V,Vi{D) (31) 
- 9jk^i = H,Vj {D) - H,V^{D). (32) 



10 



The classical Helmholtz conditions for the multiplier are recovered when we put D = 0, 
of course. 

Let us now turn to the case of forces of gyroscopic type as in (3). 

Proposition 2. The second-order field F represents a gyroscopic system of type (3) if 
and only if there is a regular function L G C°°{TQ) and a basic 2-form uj G /\^{Q) such 
that 

VOl = d^'L + ZTW. (33) 

Proof. The proof is straightforward, by a simple coordinate calculation or an argument 
like that in Proposition 1. □ 

As a preliminary remark: it is easy to verify in coordinates that for a basic 2-form uj, we 
have 

D^ZTu; = a;, d''iTUJ = 2oj, d^'iTd^uj = Sd^'oj. (34) 

It follows by taking a vertical exterior derivative of (33) that this time d^Oz will not 
vanish but must be basic, specifically we must have 

d"eL = to. (35) 

Theorem 2. The second- order field T represents a gyroscopic system of type (3) if and 
only if there exists a basic 2-form u G /\^{Q) and a symmetric type (0, 2) tensor g along 
T such that D^g is symmetric and g and co further satisfy 

Vg = 0, (36) 
^Ag-{^Agf = iTd"u, (37) 

where $ is the Jacobi endomorphism ofT. 

Proof. Assuming we are in the situation described by Proposition 2, wc define g as before 
hy g = D^'D^'L, or g{X, Y) = D^i9l(1") = D^6'l(X), from which the usual symmetry of 
D^^ follows. Acting with V on ^ and following the pattern of the proof of Theorem 1, 
we get Vg{X,Y) = d''9L{Y,X) + D''iTUj{X,Y), which is zero in view of (34) and (35). 
Finally, for the horizontal exterior derivative of (33), the modifications are that the left- 
hand side produces a term Vw in view of (35), while the second term on the right gives 
d^iTOJ = dl^co — iTd^u, and since V = on basic forms we end up with the relation 

idy^(^L - d^i^Oi = iTd^uj, 
which is the desired result (37) in view of Lemma 1. 

For the sufficiency, we observe as before that ^ is a Hessian, say g = D^D^F, and that also 
D"g will be symmetric, which in exactly the same way implies that d^Op is basic. Still 
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following the pattern of Theorem 1, = will now imply that (5 := VOp—d'^F—i'Td'^OF 
is a basic 1-form. In computing d^d^Op, the modification is that 

d^'d^ep = Id^'iav^ep = y^d^'i^ep + iTd^uj) = d"oj, 

in view of the last of (34). Since d^9p and uj arc basic, this expresses that their difference 
is closed and thus locally exact: d^Op = uj + d"a for some basic 1-form a. The compu- 
tation of d" (5 leads as before to the conclusion that j3 is closed, thus locally /3 = d" f for 
some function f on Q. Using this double information, we find that 

iTd^Op = iTOj + Va — d^iTCx, 

and subsequently 

= VOp - d"F - iTd"6p - d"f 
= V{ep -a)- d"{F - ixa + f) - itoj. 
This is a relation of type (33), with L = F — i^a + f, which concludes the proof. □ 



In coordinates, in comparison with the dissipative case of Theorem 1, the conditions (31) 
and (32) are replaced in the gyroscopic case by 

T{gij) = gikV'] + gjkV'^ (38) 



with Uij{q) = -Uji{q). 

Remark: when duj = 0, the conditions of Theorem 2 reduce to the standard Helmholtz 
conditions for a multiplier g. This should not come as a surprise, since the local exactness 
of oj then implies that the gyroscopic forces are actually of the type of the Lorentz force 
of a magnetic field, for which it is known that a generalized potential can be introduced 
to arrive at a standard Lagrangian representation. 

It is worth noting that in the sufficiency part of the proof the condition = 0, given 
that g and D^g are symmetric, is used to show the existence of a basic 1-form /3 such 
that V6p = d" F + i^d" 9p + P , where d'^Op is a basic 2-form. The condition involving $ 
then has the role of ensuring that F can be modified by the addition of a function affine 
in the fibre coordinates so as to eliminate the /3 term. This suggests that it might be 
interesting to examine the effect of ignoring the $ condition. When we do so we obtain 
the following result. 

Proposition 3. For a given second-order field T, the existence of a non-singulsr sym- 
metric type (0, 2) tensor g along r such that T)^g is symmetric and Vg = is necessary 
and sufficient for there to be a regular function L, a basic 1-form /3 and a basic 2-form 
UJ such that VOl = d"L -\- -\- j3, that is to say, such that the equations 

d f dL\ dL u ^ , , 

I - ^ = ^^ki[q)q + Pi[q), ujki = -u}ik- 



dt \dq^ ) dq 
are equivalent to those determined by T 
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Proof. It remains to show that Vg = still holds when V6l = d^L + ^t"^ + Since P 
is basic, = 0, from which it follows easily that both of the formulas cI^Ol = ^ and 
Vg{X, Y) = d"9L{Y, X) + D^ZTa;(X, Y) continue to hold, so that = as before. □ 

One point of interest about this result is that it concerns a subset of the full Helmholtz 
conditions. Unlike Theorems 1 and 2 above, but like the full Helmholtz conditions, it 
involves conditions on the multiplier only, and in this respect it anticipates the results 
to be found in the following section. 

An analogous result has been obtained by different methods in [19] (see Proposition 3.13 
and the immediately following remarks). This is the partial result that we mentioned in 
the introduction. 

4 Reduction to conditions on the multiplier only 

We have seen in the previous section that Theorems 1 and 2 produce the direct analogues 
of the Helmholtz conditions (15) of the standard inverse problem of Lagrangian mechan- 
ics. It is quite natural that the extra elements in our analysis, namely the function D, 
respectively the 2-form lo, make their appearance in these covering generalizations. Quite 
surprisingly, however, one can go a step further in the generalizations and eliminate the 
dependence on D or a; all together, to arrive at necessary and sufficient conditions in- 
volving the nuiltiplicr g only. This is what we will derive now, but it is a rather technical 
issue, for which we will therefore prepare the stage by proving a number of auxiliary 
results first. We recall that, as in the relation (9), a notation like YlxYZ what follows 
always refers to a cyclic sum over the indicated arguments. 

Lemma 2. If F e C^{TQ), Op = F and g = D^D^F then 

dlOpix, Y,z)= Yl aiRix, Y), Z), (40) 

X,Y,Z 

d%9F{X,Y,Z) = ^Rix,YMZ). (41) 

X,Y,Z 

Proof. In view of the fact that d^d^ = 0, d^Op reduces to d^iRd^F, and using (9) we 
then get 

d^j,eF{X,Y,Z)= J2 D^(zfld^F)(F,Z)= ^ (zDv^d^F + ZflD^d^F)(y,Z) 

X,Y,Z X,Y,Z 

= Y [d''F{B^xR(Y,Z))+9{X,R{Y,Z)) 

X,Y,Z 

= Y 9{RiX,Y),Z) + d''Fid''R{X,Y,Z)). 

X,Y,Z 
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Taking into account the fact that 3dXR = dXdX^ = 0, the first result follows. For the 
second there are two terms to compute. For the first we have 



X,Y,Z 



X,Y,Z X,Y,Z 



since the first term of the second line vanishes in view of the Bianchi identity d^R = 
[20]. Secondly, 

iRd"eFiX,Y,Z)= d^eFiR{X,Y),Z)= Y [^r(x,yMZ) -BpF{RiX,Y)) . 

X,Y,Z X,Y,Z 

Adding these two expressions gives the desired result (41). □ 
Lemma 3. If T)"g is symmetric then 

d«($J5-($J5f)(X,y,Z)= Y 9{yR{X,Y),Z). (42) 



Proof. We have 

d"{^jg-i^jgf){X,Y,Z)= ^ g - {^j gf){Y, Z) 

X,Y,Z 

= J2 m'^Jg + M^g-(D^<^jgf -{<PjB^gf)iY,Z) 

X,Y,Z 

= J2 [g{'Dmy),Z)-g{B^^Z),Y)\+ Y [B^gmZ)-B^g{^Z,Y) 



Making use of the cyclic sum freedom in the second and fourth term, and of the symmetry 
of D^g in the third, the right-hand side reduces to 

Y 5(Df$(y)-D?cI>(X),Z)= J2 9{d"HX,Y),Z), 

X,Y,Z X,Y,Z 

which proves our statement in view of = Vi?. □ 
Lemma 4. // g and D^g are both symmetric then 

[V, B"]g {X, Y, Z) - [V, T)"]g (F, X, Z) 

= D|($j5-($j5)^)(X,y)- Y 9iR{X,Y),Z). 

X,Y,Z 
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Proof. The commutator [V,D"] is rather comphcated when it comes to its action on a 
symmetric type (0, 2) tensor g. It reads (see for example [9] where it was already used): 

[V, D«]5 {X, Y, Z) = Bl^g{Y, Z) - 2g{R{X, Y),Z)- 2g{R{X, Z),Y) 

+ g{B^^Y),Z) + g{B^^Z),Y). 

Subtracting the same expression with X and Y interchanged, it is however a fairly simple 
computation, using the symmetry of D^g and the property = 3R, to arrive at the 
desired result. □ 

Lemma 5. For all F G C°°{TQ) we have 

D^D^D^F(X, Y, Z) - D"D^D^F(y, X, Z) = Bld"d^F{X, Y). (43) 

Proof. This is in fact a variation of a formula which was already used in proving that 

d^Op is basic in the second part of the proof of Theorem 1. We have to appeal again 
to the general formula (28), applied to D^F = d"^ F = 9p. After swapping the last two 
arguments in each term on the left in (43), a direct application of this formula easily 
leads to the result. □ 

Theorem 3. The second- order field T represents a dissipative system of type (2) if and 
only if there exists a symmetric type (0, 2) tensor g along r such that both B^g and B"g 
are symmetric and 

9{RiX,Y),Z) = Q. (44) 

X,Y,Z 

Proof. Assume first that the conditions of Theorem 1 hold true. So g and D^g are 

symmetric and as before, since V preserves the symmetry of B^g and also D^Vg' = 
D^D^D^Z) is manifestly symmetric, we conclude that D"g is symmetric. Moreover, if 
F is any function such that g = D^D^F, we know from Lemma 1 that 

^jg - ($J g)^ = idV^Op - d^i<^9F- 

It then follows from the last condition in Theorem 1 that 

= d''d''d"D = d^i^v^ep = Sd'^iOF, 

so that the first statement in Lemma 2 implies (44). 

For the converse, symmetry of g and B^g imply that Vg and VD^gi are symmetric, and 
since in addition T)"g is symmetric, we conclude that D^V^ is symmetric, which means 
that Vg is also a Hessian (see Lemma 1), say Vg = D^D^D for some function D. Next, 
we look at the statement of Lemma 4 in which the last term vanishes here by assumption. 
We have that VD^^ is symmetric, so that the left-hand side reduces to 

-D^D^D^D(X, Y, Z) + D^D^D^5(y, X, Z). 
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Combining the results of Lemma 4 and Lemma 5 we conclude that the 2-form 

is basic. Now from the last of the properties (10) and Lemma 2 applied to V^, which is 
determined by ^/j, we know that 

d^d'^d^'D = d^eD= V5(i?(X,y),Z). 

X,Y,Z 

This in turn, making use also of the result of Lemma 3, gives rise to the following 
calculation: 

d"^= E 5(Vi?(x,y),z)+ J2 V5(i?(x,y),z) 

X,Y,Z X,Y,Z 

= V( ^ 9{R{X,Y),Z)) 

X,Y,Z 

- ^ g{R{VX,Y),Z)- ^ g{R{X,VY),Z)- ^ g{R{X,Y),\7Z) 

X,Y,Z X,Y,Z X,Y,Z 

= -J2 [9iR{^Z,X),Y) + g{R{Y,VZ),X) + g{R{X,Y),VZ)]. 

X,Y,Z 

The expression between square brackets in the last line is zero because of (44), with 
vector arguments X, Y and VZ; it follows that /3 is closed, thus locally /3 = d"a for 
some basic 1-form a. Putting D = D — ixa, we have = D^D^D = D^D^D, and 
P - d"d^D = -d"d^D, so that $J 5 - ($J 5)^ = d^d"D and all conditions of Theorem 1 
are satisfied. □ 

The results of Theorem 3 deserve some further comments. Establishing necessary and 
sufficient conditions for the existence of a Lagrangian is in a way the easy part of the 
inverse problem; the hard part is the study of formal integrability of these conditions, for 
which a number of different techniques exist (see for example [2], [15, 16], [28]). If we go 
back to the standard Helmholtz conditions (15), for example, two of the first integrability 
conditions one encounters are the symmetry of D^^ and the algebraic condition (44) . So 
in the standard situation, if a 5 exists satisfying (15), these properties will automatically 
hold true: it seems to us noteworthy that these two integrability conditions make their 
appearance in the dissipative case as part of the starting set of necessary and sufficient 
conditions. It is further worth observing that the case of Rayleigh dissipation can be 
characterized by the further restriction that D^V^ = 0. Indeed, since Vg = D^D^D, 
this extra condition will imply that D must be quadratic in the velocities. 

The coordinate expressions of the conditions in Theorem 3 are, apart from (30), 

Hiigjk) - Hj (gik) + guTji^ - 5^7 = (45) 
9ijRii+9ljRik + 9kjRji = 0, (46) 
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where r/, = Vk^) and ^ = Hj{T'1) - iJ,(r^^) = - Vj{^l)). 

We now proceed in the same way for the gyroscopic case. 

Theorem 4. If the second-order field F represents a gyroscopic system of type (3) then 
there exists a symmetric type (0, 2) tensor g along r such that G^g is symmetric, Vg = 
and 

($J g - ($J gf) {X, Y)= g{R{X, Y),T). (47) 

X,Y,T 

The converse is true as well, provided we assume that $j g is smooth on the zero section 
ofTQ^Q. 

Proof. Assume we have a g and cu satisfying the conditions of Theorem 2. Acting on 
the condition (37) with d^ , the left-hand side reduces, as in the proof of the preceding 
theorem, to Sd^Op for any F such that g = D^D^F. For the right-hand side, we get 
d^iTd^cj = 3d"uj. Hence d"uj = c/^^f, and (37) can be written as 

g - i^j gf = iTd'^9F. (48) 

Making use of Lemma 2 the result now immediately follows. 

Conversely, (47) obviously implies (48) for any F such that g = D^D^F. D^g symmetric 
and = imply that D"^ is symmetric and then also [V, D"]^ is symmetric. It follows 
from Lemma 4 that 

BU^jg-{^jgf){X,Y)= ^ 5(i?(X,y),Z), yX,Y,Z. 

X,Y,Z 

In particular, taking Z to be T and using Lemma 2 again plus (48), we obtain 

D^($j g - ($j gf) = $J 5 - ($J gf- 

This asserts that ^Ag~ ($J gf is homogeneous of degree 1 in the fibre coordinates. The 
additional smoothness assumption then further implies linearity in the fibre coordinates, 
so that there exists a basic 3-form p such that <I>j — {^Agf = ztP- There are two 
conclusions we can draw from this by taking appropriate derivatives. On the one hand, 
taking the horizontal exterior derivative and using Lemma 3 we obtain 

9{yR{X, Y), Z) = {d"iTp){X, Y, Z) = (Vp - iTd"p){X, Y, Z). 

On the other, knowing that D^^tP = izP for any Z and appealing once more to the 
general conclusion of Lemma 4, wc see that actually p{X^ Y, Z) = Ylx Y z si-^i-^^ ■^)) 
from which it follows in view of = that Vp(X, Y, Z) = Y.xy,z f (vi?(^, 1"), ^). The 
conclusion from the last displayed equation is that iTd" p = 0. But then = D^ixd^p = 
ixd" p for all X, so that d" p = and locally p = d"uj for some basic uj. It follows that 

$j 5 - ($J gf = iTd"oj, 

which completes the proof. □ 
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In contrast with the preceding theorem, the condition (47) which makes its appearance 
here is not one which is directly famihar from the integrabihty analysis of the standard 
Helmholtz conditions. But indirectly, when a; = 0, the left-hand side vanishes and 
the fact that this is also the case for the right-hand side follows from the integrabihty 
condition (44). 

The coordinate expressions of the conditions in Theorem 4, in addition to (30) and (38), 
are 

SlM - 9kM = idijRii + gijRik + 9kjRiy. (49) 

Before embarking on examples, it is worth emphasizing a fundamental advantage of our 
intrinsic approach: we are not restricted to the coordinate expressions in natural bundle 
coordinates listed so far, if there are good reasons to work in a non-standard frame. This 
is the case, for example, in applications where it is appropriate to work with so-called 
quasi-velocities. Quasi-velocities are just fibre coordinates in TQ with respect to a non- 
standard frame {Xi} of vector fields on Q (which also constitute a basis for the module 
of vector fields along r). All conditions we have encountered so far may be projected 
onto such a frame and rewritten in terms of the quasi-velocities. For example, take the 
condition (19) we started from in the preceding section. It can be expressed as follows: 

= {VOl- d"L- (fD){Xi) 
= T{9l{X,)) - 9L{VXi) - Xf{L) - XUD) 
= nxUL)) - {VX.riL) - Xf{L) - X^D) 
= T{Xj{L))-Xf{L)-XnD), 

where Xf stands for the complete lift of the vector field X^. Quasi-velocities w'^ can 
be thought of as the components of T with respect to some anholonomic frame {Xi) of 
vector fields on Q. One can show (see e.g. [7]) that the complete and vertical lifts of such 
a frame, expressed in the coordinates {q,w), take the form 

where X^ = Xjd/dq^ and [Xi,Xj] = A\-Xk. The condition (19) now becomes 

^\-X^^ + A^ vj^^ = ^ 

These are the Boltzmann-Hamel equations referred to in [17], where, since the results the 
same authors obtained in [18] were expressed only in standard coordinates, all conditions 
had to be rederived from scratch. Needless to say, one can also recast any of the other 
coordinate-free conditions we have obtained in terms of quasi-velocities. 
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5 Illustrative examples and concluding remarks 



We start with a simple linear system with two degrees of freedom, which will serve us 
well to illustrate a number of features of the results we have obtained. Consider the 
system 

qi = -aqi - bq2 - wgi, (50) 
q2 = bqi - aq2 + ujq2, (51) 

where a, b and to are constant, non-zero parameters. The only non-zero connection 
coefficients are 

i 1 — 2<^ — — J- 2) 

and we obtain 

$1 = $2 ^ a - $2 = ^ = -^1- 

Since $ is constant, the curvature tensor R is zero so that condition (44) is satisfied (in 
fact it is void anyway in view of the dimension). It follows from Theorem 3 that any 
constant symmetric g should be a multiplier for a representation of the given system in 
the form (2). We consider three such non-singular matrices: 

For g^^\ it is easy to verify that with 

= U'ii - ^2) - h^(.Qi - 02) - bqiq2, 
L»i = -lio{qf + g), 

we have a representation of the given system in the form (2). In the case that a; = 0, 
5(1) is still a multiplier for the standard inverse problem and Li then becomes a genuine 
Lagrangian. Also g^'^\ which changes the order of the equations, is a multiplier in 
that case, leading to an alternative Lagrangian for the same reduced system. But that 
Lagrangian cannot serve for a dissipative representation of the full system. Instead, we 
have to take 

L2 = qm - aqiq2 - ^6(gi - qf) + ^uj{qiq2 - q2qi), 

and then D2 = 0. We discover here that the given system is variational, with L2 as 
Lagrangian. For g^^^ the situation is different again. This time, this is not a multiplier 
for the reduced system (w = 0) , it violates the condition that $j g must be symmetric. 
But for the full system, we can simply take a kinetic energy Lagrangian and then make 
a suitable adaptation for D. Explicitly, 

D3 = -a{qiqi + 92^2) + b{qiq2 - q2qi) + ^(^{qj - qf)- 
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Let us now look at the same system from the gyroscopic point of view. Since R = 
0, the rather peculiar condition (47) of Theorem 4 reduces to the usual ^-condition 
and Theorem 4 simply states the standard Helmholtz conditions for the existence of a 
multiplier. In other words, any multiplier for a representation in the form (3) will be 
a multiplier for a variational description as well. Of the non-singular, constant g^^^ we 
considered before, only g^"^^ satisfies the conditions now, and we can take 

L4 = qiq2 - aqiq2 - \h{ql - qf), 

with the 2-form co dqi A dq2 to satisfy the requirements of Theorem 2. It should of course 
not come as a surprise that we must have a variational formulation here as well, since 

we are in the situation described in the remark towards the end of Section 3. In fact we 
have already found the Lagrangian for this variational formulation: it is the function L2. 



For a second example, with n = 3, consider the non-linear system 



qi 

92 
93 



Prom (12) one easily verifies that 

(*5) = 



-2 

49293 



-9193 
^929193 + f 92^^9192 



929193, 

9i, 

9i - 9^^9293- 



-i9i93 

2^2 93 

-K''/293-i92^9t 



(52) 
(53) 
(54) 



-1^2 



and the curvature tensor R = Id^ ^ is given by 



1929193 + 19192 \ 
-^9^^9293 + 9i 
-5929? + i92"'9i J 



R = - (^Iqsdqi A dq2 - (|g2 + |9293)c^9i A dq^ - jqidq2 A dq3^ ® ^ 
+ (591^^91 A dqs - \q2^qidq2 A dgs^ (g) — 
- (19^^91 c^9i A dq2 + \q2q1dq1 A ^^3 - |9^^92C?92 A ^^3^ (g) 



The multiplier problem is already quite complicated for a system of this kind and it is 
not our intention here to explore all possible solutions. For simplicity, therefore, we limit 
ourselves in the dissipative case (2) to analyzing the existence of a diagonal multiplier g 
which depends on the coordinates qi only. With such an ansatz, the curvature condition 
(44) in Theorem 3 reduces to 

933 = (911 - 2922)92, 

and the requirement that g should be symmetric subsequently imposes that gu = 
4^22 = constant. Hence, up to a constant factor, we are reduced to the possibility that 



511 



522 



1, 



533 = 2^2 • 
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As was mentioned in the previous section, the conditions imposed so far are also integra- 
bihty conditions in the standard inverse problem so that, starting from the same ansatz, 
this g would also be the only candidate for a standard Lagrangian representation of the 
system. When we compute Vg now, we get 



((V5)«,) 



4^293 4^291 


4^291 



Since Vg ^ 0, our candidate cannot lead to a variational formulation. On the other 

hand. Theorem 3 is satisfied, so there must exist a D for a dissipative representation. 
From the requirement (24) in Theorem 1, one easily verifies that such a D must satisfy 

= 4^29193 + /ll, y2{D) = h2, V3{D) = 2q2ql + h, 

where the hi are as yet arbitrary functions of the coordinates. The final requirement 
(25) of Theorem 1 then shows that the hi can be taken to be zero. Thus, 

L=^(4ql+g + 2q2ql) and D = 2q2qfq3, 

provide a solution for the inverse problem of type (2) for the given system. 

Concerning the inverse problem of type (3), it is less appropriate to look for a diagonal 
g (as the example with n = 2 has shown), but even if we extend our search to a general 
g depending on the qi only, the conditions of Theorem 4 have no non-singular solution. 

Consider, finally, the system 

qi = bgi94, (55) 

92 = 9294, (56) 

93 = (1 - b)qiq2 + 69291 94 - ?'9i9294 + {b+ l)g394, (57) 

94 = 0, (58) 

with —l<b<l and b ^ 0. These equations can be interpreted as the geodesic 
equations of the canonical connection associated with a certain Lie group G, which is 
uniquely defined by VxY = ^[-^, where X and Y are left-invariant vector fields. In 
the case of the above system, the Lie group is listed as ^4^9^ in [12], and it was shown (see 
also [1]) that the system does not have a variational formulation. This is a consequence 
of the integrability condition (44) which can only be satisfied by multipliers for which 
913 = 923 = 533 = 0. But then, the <I>-condition in (15) leads automatically to 534 = so 
that there is no non-singular solution. 

Notice that the system is invariant for translations in the q^ and 54 direction; it is there- 
fore reasonable that we limit ourselves in our search for non-conservative representations 
to multipliers with the same symmetry. In the dissipative case, after haven taken the 
same curvature condition into account, the D^^-condition leads to the further restrictions 

^^34 „ 9^34 „ 

-^T- = 0, 534 = 0, 

dqs dqz 



21 



among others, from which again 334 = follows, with the same negative conclusion. 
In the gyroscopic case, one can show that the condition (47) cannot be satisfied for a 
multiplier with coefficients depending on the coordinates qi only. 

Some final comments are in order. In the case of linear systems such as our first example, 
it frequently happens that a multiplier for the inverse problem exists which is a function 
of time only (see for example [25] ) , so we will briefly sketch here how our present theory 
can be extended to general, potentially time-dependent second-order systems. First of 
all, the extension of the calculus along r : TQ — ?• Q to a time-dependent setting has 
been fully developed in [29], which also contains the analogues of the conditions (14) 
and (15) for the inverse problem. In all generality, we are then talking about a calculus 
of forms along the projection tt : M X TQ M x Q say. But as has been observed for 
example in [9], the extra time-component in this setting docs not really play a role when 
it comes to studying the Helmholtz conditions and their integrability. That is to say: 
one has to use dt and the contact forms dg* — q^dt as local basis for forms along tt and 
a suitable dual basis for the vector fields which includes the given second-order system 
T; important geometrical objects such as the Jacobi endomorphism $ will pick up an 
extra term for sure, but when restricted to act on vector fields without F-component, 
all formulas of interest formally look the same. It is therefore not so hard to apply a 
suitably reformulated version of the present theory to time- dependent systems, when 
needed. That important formulas formally look the same will be seen also in the final 
observations in the appendix, where the setting is essentially time-dependent, though the 
approach adopted there is quite different again from the calculus along tt we are referring 
to here. 



Appendix 

In this appendix we will relate our results to those obtained by Kielau et al. in [17, 18], 
especially the latter; but first we wish to derive those results anew, in a way which allows 
us to explain an interesting feature of them which was mentioned in [18] but not fully 
dealt with there. 

The problem discussed in [17, 18] differs in several ways from the one which has been the 
subject of our paper, the most important of which is that it is assumed there that a system 
of second-order ordinary differential equations is given in implicit form fi{t,q,q,q) = 0, 
and the problem posed is to find necessary and sufficient conditions on the functions fi 
such that the equations may be written in the form (2), where L and D are allowed to 
be time-dependent. That is, the question is whether the equations are of Lagrangian 
type with dissipation as they stand, rather then whether they may be made equivalent 
to such equations by a choice of multiplier. 

It is probably most satisfactory to approach the inverse problem for a second-order system 
given in implicit form by using the methods associated with variational sequences, rather 
than the techniques employed in the body of the paper. Fortunately we will need only 
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the rudiments of such methods, one version of which we now briefly describe; justification 
for the unsupported claims we make can be found in [31], for example. 

We deal with the (trivial) fibred manifold tt : Q x M — )■ M, and its infinite jet bundle 
J°°(7r); this may seem a bit extravagant when we are interested only in second-order 
equations, but is convenient for technical reasons. However, all functions and forms 
under consideration will be of finite type (i.e. depend on finitely many variables). We 
take coordinates (t,^*) on R x Q; the jet coordinates are written g*, and so on. We 
denote the contact 1-forms by 

9' = dq' - q'dt, 9' = dq' - q'dt, 9' = dif - q'dt, .... 

We need two exterior-derivative-like operators on exterior forms on J°°(7r). The first is 
the vertical differential dy (not to be confused with d^), which is defined by 

= + ^i^' + 14^' + • • • ' ^ydt = dvO' = dv9' = dv9' = . . . = 0. 
ag* oq^ oq^ 

The key properties of dy are that dy = 0, and that dy is locally exact. The second 
operator is the variational differential (5, about which we need to say just the following. 
First, for a Langrangian L(i, q, q) 

5{Ldt) =Ei{L)9' Adt 

where the Ei[L) arc the Eulcr-Lagrangc expressions. The 2-forms which, like 5{Ldt), are 
linear combinations of the 0* A dt are called source forms in [31] and dynamical forms in 
[19]. Secondly, for any source form e = A dt, fi = fi{t, q, q, q), 

Se = -\{rij9' A 9^ + Sij9' A 9^ + tij9' A 9^) A dt, 

where the coefficients are given by 

i±(dh_df£\ .^(dh_dfj\ 
2 dt \dq^ 2 ^^2 yQqj Qqi ] 

dt \dqi J 



Again, 6^ = and 6 is locally exact. With a source form e = fi9'^ A dt, fi = fi{t, q, q, q), 
one associates the second-order system /j = 0, and conversely; so that (5e = is necessary 
and sufficient for the second-order system /j = to be locally of Euler-Lagrange type. 
The vanishing of the coefficients rij, Sij and tij are the (classical) Helmholtz conditions 
(see e.g. [18, 19, 24]). 

By considering the transformation properties of the jet coordinates and the contact forms 
under transformations of the form q^ = q^{t,q), t = t one can show that the set of forms 



dqi dq' 

dq^j dq^ 
dip dq^ 
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spanned by {dt,9^,0^} with coefficients which are functions of t, q and q is well-defined. 
We call forms like this first-order forms. Note that dy maps first-order forms to first- 
order forms. Moreover, one proves exactness of dy by using essentially the homotopy 
operator for the de Rham complex for the variables q, q, . . . , treating t as a parameter. 
It follows that if a is of first order and satisfies dyct = then there is a first- order form 
/3 such that a = dyP- 

The first step in applying these concepts to dissipative systems is to characterize dissi- 
pative force terms using them. 

Proposition 4. The first-order source form 

A = A dt 
dq' 

satisfies SA = dyA. Conversely, if £ is a first-order source form such that 5e = dye 
then £ = A for some first- order function D. 

Proof. Applying the formula for 6 acting on a source form one finds that 

6A = -^^——Q' AO^ Adt- ^-^9' AO^ Adt 
oqWq'- oq^oq^ 

= dyA. 

Conversely, if £ = fiO^ A dt is of first order then tij = and 

dq^ dq^ 

But 

dye = ^9' A 9^ Adt- ^9' A 9^ A dt, 
oq^ oqO 



and so if Se = dye then 

Sin = 7TT7 + TTij = ^TTTl OV 



'•^ dq^ dq^ dqi dq^ dq^ ' 



so that there is a function D = D{t, q, q) such that fi = dD/dq^. The terms in 9^ A9^ Adt 
then agree. □ 

Now take e to be the source form representing the given equations. If they arc of Euler- 
Lagrange type with a dissipative term then e = 5{Ldt) — A, so 5e = 5 A = dyA. Then 
Se is of first order, and furthermore dy6e = 0. These are necessary conditions for the 
given system to take the desired form. They are in fact sufficient also, as we now show. 
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Theorem 5. A system of second-order ordinary differential equations fi{t,q,q,q) = 
may be written locally as 

d_ /dL\ dL _ dP 
dt \dq^ J 5g* 

for some first-order functions L and D if and only if the corresponding source form e is 
such that 5e is of first order and satisfies dySe = 0. 

Proof. It remains to prove sufficiency. By the local exactness of dy we may assume 
that Se = dya for some first-order 2-form a (not necessarily a source form), which 
is determined only up to the addition of a dy-exact form. Contact 2-forms fi can be 
ignored in a since dyfi then contains no dt terms and therefore cannot contribute to Se. 
Moreover, if we put 

a = {Xiff'- + jj.iff'-) A dt, 
since de contains no 6^ A 9^ terms either, we must have 

dq^ dq^ ' 

so that Hi = dil)/dq^ for some function ^(t, g, q). Hence, up to the dy-exact term dy{'tpdt), 
a is of the form 

a = - A dt, 

i.e. is a source form, say a = Ui6^ A dt. Since Se is of first order we must have tij = 0, 
whence Sij is symmetric in i and j. It follows that the coefficient of 0* A 6^ in dya 
must be symmetric in i and j, whence Vi = —dD/dq^ for some function D = D{t,q,q), 
and a = —A. Then Se = dya = —dyA = —SA, and so there is some L such that 
£ = S{Ldt) — A, as required. □ 

The point of interest that we mentioned at the beginning of this appendix is that this 
version of the result represents the conditions in part as the closure (under dy) of a 
certain form, namely Se. That it might be possible to state the conditions in such a way 
was raised speculatively in [18], but the form and operator were not specifically identified 
there. 

Under the assumption that Se is of first order, so that tij = 0, it is easy to verify that 
the o?y-closure conditions are 



= 




dq'- 




dq^ 


dqi 


drij 


dsik 


dsjk 




dq^ 


dqi 


dq^ 




dsij _ 


dsik 






dq^ 


dqi ' 
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It must not be forgotten that r^j and Sij are supposed to be of first order; in addition, 
tij = means that we have 

dip dq^ 

These are the generahzed Helmholtz conditions as given in [18]. However, it turns out 
that the first and last of the closure conditions are consequences of the other conditions, 

as wc showed in [8]. It then follows easily that the following conditions are equivalent to 
those given above: /j = Qijq^ + hi with gij symmetric, where gij, hi are of first order and 
further satisfy 





dgik 




dq^ 


dq^ 




dgik 


1 d'^hi _ dgjk 


1 9% 


dqi 


^ dq^dq^ dq^ 


^ dq^dq^ 


El 

A A I, 


f d'^hi d'^hi \ 


= 0, 


\dq^dq^ dq^dq^ J 



where ^ ^ stands for the cyclic sum over the indices. 

As the problem has been presented so far in this appendix, we must take gij and hi as 
given; the equations above provide a test for determining whether the given second-order 
system can be put into the required form. However, it is now possible to regard these 
equations from the alternative point of view: we set hi = gijf^ , where the second-order 
system is given in normal form q^ = P{t,q,q); we regard the /* as known but the gij 
as to be determined; the equations above now become partial differential equations for 
the unknowns gij. We leave it to the reader to verify that they correspond (take the 
autonomous case for simplicity), in the order written, to the conditions of Theorem 3, 
namely D^y is symmetric, D^^ is symmetric, and y z di^i-^^^)^ ^) ~ ^■ 
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